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, Abstract. We discuss the effects of rotation on the evolution of accreting carbon-oxygen white dwarfs, with the emphasis 

' on possible consequences in Type la supernova (SN la) progenitors. Starting with a slowly rotating white dwarf, we consider 

' the accretion of matter and angular momentum from a quasi-Keplerian accretion disk. Numerical simulations with initial white 

, dwarf masses of 0.8, 0.9 and 1.0 Mq and accretion of carbon-oxygen rich matter at rates of 3 ... 10 x 10~^ Mq/jv are per- 

fS) ' formed. The models are evolved either up to a ratio of rotational to potential energy of T/W — 0.18 — as angular momentum 

loss through gravitational wave radiation will become important for T/W < 0.18 — or to central carbon ignition. The role of 
the various rotationally induced hydrodynamic instabilities for the transport of angular momentum inside the white dwarf is in- 
vestigated. We find that the dynamical shear instability is the most important one in the highly degenerate core, while Eddington 
Sweet circulation, Goldreich-Schubert-Fricke instability and secular shear instability are most relevant in the non-degenerate 
envelope. Our results imply that accreting white dwarfs rotate differentially throughout, with a shear rate close to the threshold 
value for the onset of the dynamical shear instability. As the latter depends on the temperature of the white dwarf, the thermal 
, evolution of the white dwarf core is found to be relevant for the angular momentum redistribution. As found previously, signif- 

' leant rotation is shown to lead to carbon ignition masses well above 1.4 Mq. Our models suggest a wide range of white dwarf 

I explosion masses, which could be responsible for some aspects of the diversity observed in SNe la. We analyze the potential 

. role of the bar-mode and the r-mode instability in rapidly rotating white dwarfs, which may impose angular momentum loss 

by gravitational wave radiation. We discuss the consequences of the resulting spin-down for the fate of the white dwarf, and 
the possibility to detect the emitted gravitational waves at frequencies of 0.1 . . . 1.0 Hz in nearby galaxies with LISA. Possible 
implications of fast and differentially rotating white dwarf cores for the flame propagation in exploding white dwarfs are also 
O ' briefly discussed. 

^ ' Key words, stars: white dwarf - stars: accretion - stars: rotation - supernovae: Type la - gravitational waves 
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1 . Introduction the Phillips' relation to very distant SNe la. An understanding 

5J] of the origin of the diversity observed in SNe la is thus a crucial 

. Type la Supernovae (SNe la) have a particular importance in ^^^^ ^^^jj^. evolution theory, which requires to identify the 

astrophysics. Observations of SNe la at low redshift showed ^^^^^^^ evolutionary paths of SNe la progenitors, 
a clear coiTelation between the peak brightness and the width 

of the hght curve (PhilUps i993. Phillips' relation), which al- Unlike core collapse supernovae. Type la supernovae (SNe 

lowed to use SNe la as distance indicators for galaxies even ^""e believed to occur exclusively in binary systems (e.g. 

beyond z = 1. This made SNe la an indispensable tool for cos- ^ivio 2001 ). Although it is still unclear which kinds of bi- 

mology, in particular to determine the cosmological parameters ^y^^ems lead to SNe la, non-degenerate stars such as main 

(e.g. Hamuy et al.|l996l Branch[T998| Leibundgut|2000||200D. sequence stars, red giants or helium stars are often assumed 

The new cosmology with a non-zero cosmological constant has white dwarf companion (e.g. Li & van den Heuvel 

been initiated by the observational evidence deduced from ob- "^'^hisu et al. 1999 ; Langer et al. 2000 Yoon & Langer 2003). 

servations of SNe la at high redshift (Perlmutter et al. ^^'^^ "^^^'^ ^'''''^^'^ ^° ^° Chandrasekhar 

Riess et al I SOOO t Umit by mass accretion from its companion, with accretion 



, j-r.xTTi 111 if^xT rates which allow steady shell hydrogen and helium burning 

Recent analyses of SNe la have revealed, however, that SNe ,,V\i^-.7i,r / n» 1 i. ,-i i- 

. J. , , ,1 .... (Al > 10 ' JVl(7)/yr). An understanding of the physics of mass 

la are not perfectly homogeneous but show some diversity in . , , . . , , . 

, . ,. , J / n 1 i>^7wn NT accretion IS therefore indispensable to investigate the evolution 

their light curves and spectra (e.g. Branch 120011 Nomoto ^t ' h' d f 

aL l2003i. Li et al. 20031) . This leaves concerns about applying ° 

Although the mass accretion process in white dwarfs has 

Send ojfprint requests to: S.-C. Yoon (email: been discussed by many authors (e.g. Iben'1982 i Nomoto' 1982; 

s.c.Yoon@astro.uu.nl) Fujimoto & Sugimoto .1982. Saio & Nomoto .1985. .1998. 
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Kawai et al. fMS' Cassisi et al. [T9981 Piersanti et al. l2000l 
Langer et al. 2002 1, little attention was so far devoted to the ef- 
fects of angular momentum accretion and the ensuing white 
dwarf rotation (see Sect. 13. 11 1. As the evolution of stars can 
generally be affected by rotation (e.g. Heger & Langer 120001 
Maeder & Me vnet 120001 . this might be particularly so in ac- 
creting white dwarfs: Since the transfered matter from the 
white dwarf companions may form a Keplerian disk which 
carries a large amount of angular momentum, the resultant 
accretion of angular momentum will lead to the spin-up of 
the white dwarf (e.g. Durisen f 19771 Ritter 1985 Narayan & 
Popham 1989 Langer et al. '20^' 2002 l20D?t . The obser- 
vation that white dwarfs in cataclysmic variables rotate much 
faster than isolated ones (Sion 199$') provides evidence for ac- 
creting white dwarfs indeed being spun up. Rapidly rotating 
progenitors may also lead to aspherical explosions, which may 
give rise to the observed polarization of SNe la (Howell et 
al. l200n Wang et al. l2003l . 

Here, we make an effort to investigate in detail the possibil- 
ity of angular momentum accretion, and the role of the various 
rotationally induced hydrodynamic instabilities in transporting 
angular momentum into the white dwarf core, and in establish- 
ing the the pre-explosion angular momentum profile. The re- 
mainder of this paper is organized as follows. We evaluate the 
possible mechanisms for angular momentum transport in ac- 
creting white dwarfs in Sect. |2] Our approach to the problem, 
including the numerical methods and physical assumptions, is 
discussed in Sect. |3l where previous papers on rotating white 
dwarf models are also reviewed (Sect. l3JT l. Numerical results 
are presented in Sect. with the emphasis on the process of 
angular momentum transport in the white dwarf interior. Pre- 
explosion conditions of accreting white dwarfs and their impli- 
cations for the diversity of SNe la are discussed in Sect. |5] |6l 
and The possibility of detecting gravitational waves from 
SNe la progenitors is examined in Sect.ls] Our conclusions are 
summarized in Sect.|9] 



2. Angular momentum transport In white dwarfs 

Inside a white dwarf, angular momentum can be transported 
by Eddington Sweet circulations and by turbulent diffusion in- 
duced by hydrodynamic instabilities. Our discussion below is 
limited to angular momentum transport in the vertical direction. 



2. 1. Eddington-Sweet circulation 

Eddington-Sweet circulations are induced by the thermal im- 
balance between the equator and the poles of a rotating star. 
Its time scale is roughly given by tes — tkh/x^ (Maeder & 
Meynet'2000), where tkh is the Kelvin-Helmholtz time scale 
and X is the angular velocity normalized to the Keplerian value, 
i.e., X — ^I'^Kc-p- In a white dwarf which accretes at rates 
M> 10^^ Mo/yr, the surface luminosity reaches lO"*^ L0 due 
to compressional heating and nuclear burning. The quantity % 
is close to 1 in the outer envelope as we shall see in 14.3.51 As a 
result, we expect an Eddington Sweet circulation time scale in 
the non-degenerate envelope which is shorter or comparable to 



the accretion time scale: 

Tes 



RLx^ 
1.3 • lO^yr 



(AMe„v/0.05MQ)2 



(i?/O.OlRQ)(i/lO4L0)(x/O.8)2- 



(1) 



In the degenerate core, on the other hand, we have tes ^ 
Tacc since the Kelvin-Helmholtz time scale is typically much 
larger than 10^ yr. Furthermore, x is usually significantly 
smaller than 1 in the inner core: 



GMl 



Tes 



-^corc-^coreX 

7 • lOVr 



(Mcorc/Mo)2 



(i?,orc/O.OO5RQ)(L,orc/lOL0)(x/O.3)2 



(2) 

Consequently, Eddington-Sweet circulations are only impor- 
tant in the white dwarf envelope. 

2.2. Siiear instability 

The dynamical shear instability (DSI) occurs when the energy 
of the shear motion dominates over the buoyancy potential. In 
a rotating flow, the linear condition for instability is given by 



< Ri.c « 1/4, 



(3) 



where Ri c denotes the critical Richardson number which is 
about 0.25, N'^ the Brunt-Vaisala frequency, and a the shear 
factor (Zahn [T974l : 



9l_ 

Hp 



and 



(Vad-V), 7V2^ 



du! 
d\nr 



Hp 



(4) 



(5) 



(6) 



Here, u! is the angular velocity and other symbols have their 
usual meaning. The equation of state of electron degenerate 
matter is given by P = T, ^e,Mi) ™d the factor V'V^ 
in Eq. |5] needs to be replaced by (yScV^^ + (piV^tj, where 

ipc = (31n/9/91n^o)p,T,Mi and ipi = {d\n p/d\n fii)p^T,fi^- 
In a typical CO white dwarf, N'^ is generally smaller than 
and has minor effects on the angular momentum transport as 
discussed in Sect. 14.3.4) 

Under constant temperature and homogeneous chemical 
composition, which is a good approximation for the central re- 
gion of a CO white dwarf, the critical value of a for the DSI 
can be given as 



''DSI.crit - N /Ri 

9 



0.2 



109cm/s2 J VO.Ol / V8 • low 



Vad 

0.4 



(7) 



where Ri ^ — 0.25 is assumed. Fig.[2shows |(JDSi.crit| as func- 
tion of density, with other physical assumptions as indicated 
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Fig. 1. Threshold values of the shear factor (T = {duj/d\nr)foi 
the dynamical (solid line) and the secular (dashed line) shear 
instability. Constant gravity and temperature, i.e., g = 10^ 
cm/s^, T — 5 X 10^ K, are assumed. The chemical compo- 
sition is also assumed to be constant, with Xc = 0.43 and 
Xq = 0.54. For the critical Richardson number, i?i c = 1/4 is 
employed. i?e,c = 2500 is used for the critical Reynolds num- 
ber. 

in the figure caption. This figure shows that the condition for 
the DSI becomes more relaxed with increasing density, which 
is due to the buoyancy force becoming weaker with stronger 
degeneracy. This implies that the highly degenerate core in a 
white dwarf is more susceptible to the dynamical shear insta- 
bility than the outer envelope. 

The stability criterion for the shear instability can be re- 
laxed if thermal diffusion reduces the buoyancy force such that 



^is.l • — 77 9 *^ ^i,C5 



(8) 



where i?o,c is the critical Reynolds number, i/ is the kinetic 
viscosity and K — {AacT^) / {3CpKp^) is the thermal diffu- 
sivity (e.g. Maeder & Meynet 2000). Since the /i gradient is 
not affected by the thermal diffusion, this condition needs to be 
supplemented by the following additional condition (Endal & 
Sofia[T976j: 

Ris,2 ^ < Ri,c (9) 



The thermally induced shear instability is often called 'secular 
shear instability' (SSI). 

This instability can occur only when the thermal diffusion 
time scale is shorter than the turbulent viscous time scale (i.e., 
z/i?c,c < K), which is often the case in non-degenerate stars. 
However, this condition is not always fulfilled in white dwarfs. 
In Fig.|2j the critical density pssi.crit above which the secular 
shear instability is not allowed is plotted as function of temper- 
ature. In this calculation, the ion and electron viscosity as well 
as the radiative and conductive opacities are estimated as de- 




Fig. 2. The critical density pssi.crit, above which the secular 
shear instability can not occur, as function of temperature (solid 
line). The region where the thermal diffusion time scale ex- 
ceeds the turbulent viscous time scale is hatched. Abundances 
of Xq = 0.43 and Xq = 0.54 are assumed. i?c.c = 2500 is 
used for the critical Reynolds number 



scribed in Sect. 13.21 This figure indicates that the SSI may play 
a role only relatively weakly degenerate regions. 

The threshold value of the shear factor for this instability 
is given by aigi^^Ht = yR^,cN^/{KRi^^) if iV^ = q, and plot- 
ted in Fig [2 with the assumptions as indicated in the figure 
caption. As shown in the figure and as discussed above, ther- 
mal diffusion can not weaken the restoring buoyancy force for 
p ?^ 5 X 10^ g/cvo? under the chosen conditions, and only the 
dynamical shear instability can occur for higher densities. 

Zahn il992j derives the diffusion coefficient for the SSI as 



SI 



1 Ka^Ri^c 
3 iV2 



(10) 



The time scale of the SSI when a = CDSi.crit thus be- 
comes simply the thermal diffusion time scale: rssi (at cr — 
CDSi.crit) ~ 3R^/K. In the degenerate core of a CO white 
dwarf with p < pssi,crit, we have tssi « 7.6 x 10'^(i?corc/2 • 
10^ cm)V(ii:/500 cm^s^i) yr 

As consequence of these considerations, any shear mo- 
tion with (7 > CTDSi.crit will dccay such that a approaches 
CDSi.crit on the dynamical time scale. Further angular momen- 
tum transport will operate on the thermal diffusion time scale 
until <T reaches crssLcrit for p ^ pssi,crit- For a < crssi,crit 
and a < (TDSi.ciit (hatched region in Fig.[3, angular momen- 
tum will be transported only via the electron viscosity or via 
Eddington Sweet circulations on much longer time scales (un- 
less other kinds of instabilities are invoked). Therefore, if we 
consider accretion at rates required by the single degenerate SN 
la progenitor scenario (M > 10^^ M© /yr), the degree of shear 
may not be far from the threshold value for the dynamical shear 
instability throughout the degenerate white dwarf interior. This 
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conjecture is confirmed by the numerical results presented in 
Sect.lO 

2.3. GSF and magnetic instabilities 

The Goldreich, Schubert and Fricke instability (GSF insta- 
bility) can be induced if a star is in a baroclinic condition 
(Goldreich & Schubert 119671 Fricke I1968> . In an accreting 
white dwarf, this instability may be important in the non- 
degenerate envelope, but is likely suppressed in the degen- 
erate core where the barotropic condition will be retained 
through a dynamical meridional circulation (Kippenhahn & 
Mollenhoff ll974l see also Durisen'1977'l. Magnetic instabil- 
ities such as the Tayler instability (Spruit 2002 1 may be poten- 
tially important, and their role will be investigated in the near 
future. In the present study, we restrict our discussions to non- 
magnetized white dwarfs. 

3. Simulations of accreting white dwarfs 

For the detailed investigation of the angular momentum trans- 
port and its role in the evolution of accreting white dwarfs, we 
have simulated accreting CO white dwarfs with various initial 
masses and accretion rates. 

3. 1. Previous work on rotating wiiite dwarfs 

In order to compare our results with previously obtained 
ones, we briefly discuss the models of rotating white dwarfs 
in the literature here. Most rotating white dwarf equilibrium 
models have been constructed assuming a barotropic equa- 
tion of state and axial symmetry for both, rigidly and dif- 
ferentially rotating cases (e.g. James '1964; Roxburgh U965I 
Monaghan J 966. Ostriker J 966. Anand .1968. Ostriker & 
Mark [T968l Hachisu [l986> . These models have mainly been 
used to investigate the stability of rapidly rotating zero tem- 
perature white dwarfs (e.g. James 1964. Lynden-Bell & 
Ostriker [T967l Ostriker & Bodenheimer [T9^ Durisen [T975bl 
Durisen & Imammura 19811. One of the main results of these 
studies is that the Chandrasekhar limit increases slightly in 
rigidly rotating white dwarfs (^ 1.48 Mq), while differen- 
tially rotating white dwarfs can be dynamically stable even 
with masses up to ^ 4.0 Mq. It was also found that differen- 
tially rotating white dwarfs become secularly unstable to grav- 
itational wave radiation if the ratio of the rotational energy to 
the gravitational energy exceeds a certain limit (see Sect.|5lfor 
more detailed discussions). 

The evolution of rotating white dwarfs was investigated by 
Durisen il973a 1973b 1975a i for the first time. Evolution and 
spin-up of mass accreting white dwarfs were probed by Durisen 
J 197 7"). In his studies based on two-dimensional barotropic dif- 
ferentially rotating white dwarf models, Durisen assumed that 
angular momentum in the white dwarf interior is only trans- 
ported by electron viscosity, and angular momentum transport 
time scales of the order of 10^*' yr have been obtained. In all 
these models, the thermal history and the energy transport in 
the white dwarf matter were neglected, assuming zero temper- 
ature. 



Recently, Piersanti et al. J2003> investigated the effect 
of rotation on the thermal evolution of CO-accreting white 
dwarfs, using a one dimensional stellar evolution code. They 
also discussed angular momentum loss by gravitational wave 
radiation (cf. Sect. |5}- However, the detailed history of the 
angular momentum transport inside the white dwarf was ne- 
glected, but rigid body rotation was assumed. I.e., while 
Durisen assumed much slower angular momentum redistri- 
bution by restricting his considerations to electron diffusion, 
Piersanti et al. assume instant angular momentum redistribu- 
tion with a maximum efficiency. 

Uenishi et al. J2003t did consider a finite spin-up time of 
accreting white dwarfs by dividing their two-dimensional white 
dwarf models into a non-rotating central core and a fast rotating 
outer envelope. However, rigid body rotation is again assumed 
for the spun-up region, introducing a discontinuity in the angu- 
lar velocity profile. 

In the present study on the evolution of accreting white 
dwarfs, we consider the effects of the accretion induced heating 
and energy transport, the angular momentum transport by var- 
ious instabilities, and the effect of rotation on the white dwarf 
structure. However, our numerical models are calculated in a 
one dimensional approximation as described in the next sec- 
tion, and thus the structure of rapidly rotating white dwarfs is 
described less accurate than in multi-dimensional models men- 
tioned above. Furthermore, by assuming CO-accretion we do 
not consider effects of nuclear shell burning. These, and the 
feedback between rotation and nuclear shell burning, as well as 
rotationally induced chemical mixing, are discussed in a sepa- 
rate paper (Yoon et al. 2004 1. The effects of rotation with re- 
spect to the off-center helium detonation models - i.e., helium 
accreting CO white dwarf models with M =r-^ 10~^ Mq /yr - 
are investigated in Yoon & Langer (.2004i) . 

3.2. Numerical mettiod 

We use a hydrodynamic stellar evolution code (Langer et 
al. I1988t for the accreting white dwarf model calculations. 
Radiative opacities are taken from Iglesias & Rogers ( I1996> . 
For the electron conductive opacity, we follow Hubbard & 
Lampe J1969> in the non-relativistic case, and Canuto J 19701 1 in 
the relativistic case. The accreted matter is assumed to have the 
same entropy as that of the surface of the accreting white dwarf, 
and the accretion induced compressional heating is treated as 
in Neo et al ( I1977t . 

The effect of rotation on the stellar structure is approxi- 
mated in one dimension by introducing the factors /p and /t 
in the momentum and energy conservation equations (cf. Heger 
et al. 2000), following the method of Kippenhahn & Thomas 
C.1970 ) and Endal & Sofia ( 1976). This method is well suited 
for the case where the effective gravity can be derived from an 
effective potential. This is indeed a good assumption for a CO 
white dwarf, since it should be in a barotropic condition except 
for the non-degenerate outer envelope. Note that this method is 
also applicable in the case of shellular rotation as discussed by 
Meynet&Maeder(1997). 
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However, our method of computing the effective gravita- 
tional potential in a rotating star limits the accuracy of our re- 
sults for very rapid rotation. The potential is expanded in terms 
of spherical harmonics, of which we only consider terms up to 
the second order (cf., Kippenhahn & Thomas 1970 1. Fliegner 
J1993t showed this method to accurately reproduce the shapes 
of rigidly rotating polytropes up to a rotation rate of about 
60% critical, corresponding to correction factors of /p ~ 0.75 
and /t — 0.95 in the stellar structure equations (cf. Heger et 



al . 120001 1. We therefore limit these factors to the quoted values, 
with the consequence that we underestimate the effect of the 
centrifugal force in layers which rotate faster than about 60% 
critical. As in our models the layers near the surface are always 
close to critical rotation, the stellar radius of our models may 
be underestimated. Furthermore, our models are per definition 
rotationally symmetric. Therefore, we are in principle unable to 
investigate the onset of triaxial instabilities, which may affect 
the final fate of a rapidly rotating white dwarf (Sect.|5}. 

The angular momentum transport induced by the insta- 
bilities mentioned in Sect. |2l is described as a diffusion pro- 
cess (Heger et al. 2000 1, while this approximation neglects the 
possibility of advective angular momentum redistribution by 
Eddingtion Sweet circulations (Maeder & Meynet 20001. Note 
also that this approach is based on the assumption of shellular 
rotation (see Heger et al. l2000l for a detailed discussion), which 
might not be appropriate for white dwarfs. Nevertheless, our 
models can represent the case of cylindrical rotation to some 
degree, since most of the total angular momentum is confined 
to layers near the equatorial plane in both cases. Since the dy- 
namical shear instability is important in the present study, the 
diffusion solver has been improved as a non-linear process as 
explained in Yoon ( 2004 1, in order to properly deal with such a 
fast angular momentum redistribution as on a dynamical time 
scale during the secular evolution of accreting white dwarfs. 
Diffusion coefficients for each of the instabilities are taken 
from Heger et al. ( 2000) and Heger and Langer J2000t 
some modifications as follows. 



with 



As mentioned in Sect. 12.31 the GSF instability may be sup- 
pressed in the degenerate barotropic core, and we describe this 
effect as 



Dgsf = ^'gsf ( 1 - 



1, 



Pi 



ideal 



(11) 



Here, i'Qgp is the diffusion coefficient for the GSF instability 
as estimated in Heger et al. (2000). Pc is the electron pressure 
for complete non-relativistic degeneracy, and Pidcai is the ideal 
gas pressure: 



Pc = 9.91 • 10^^ 



5/3 



Pi, 



deal 



pkT 
/imn 



(12) 



For the secular instability, Zahn's J 19921 1 prescription is em- 
ployed (Eq. I10> . with a correction factor as in Heger et al. 
(2()()()i to ensure a smooth turn-on of the instability and to re- 
late its strength to the deviation from the Richardson criterion: 



-Dssi — - 



1 Ka^Ri , 



7V2 



1 - 



max{Pis,i, Pis, 2} 

Pi,c 



(13) 



Table 1. Physical quantities of the initial white dwarf models: 
mass, surface luminosity, central temperature, central density, 
radius and rotation velocity 



AfwD.init 




Tc,init 


Pc.init 


fiwD.init 


^rot ,init 


Mq 


Lq 


10® K 


10^ g/cm=' 


Rq 


km/s 


0.8 


0.80 


0.79 


0.97 


0.0115 


10 


0.9 


0.88 


0.85 


1.68 


0.0100 


10 


1.0 


0.81 


0.94 


3.14 


0.0086 


10 



where Pis,i and Pis. 2 are defined in Eq. |S]and|5] As already 
discussed in Sect. 12.21 the viscosity is an important factor for 
the stability of the secular shear instability. The former way of 
computing the ion viscosity in Heger et al. ( 2000) and Heger 



& Langer ( 2000 1, following Spitzer J 19621 1. is not well adapted 
for matter in the ion liquid state: here we follow the work of 
Wallenborn and Bauss Ll978j . as described in Itoh et al. (I1987> : 



5.53- 103^^-1/3(^/10^)5/^* (gcm-^s-i). (14) 



We constructed a fitting formula for the dimensionless param- 
eter rj* from the table in Itoh et al. as 

ri* =- 0.016321227 + 1.0198850 r-i-^^mro 

+ 0.024113535 r°-49999098 ^ (15) 

where T = (Ze)^/[(3/47rni)i/3fcr]. The electron viscosity, 
which is dominant over the ion viscosity in degenerate matter, 
is estimated according to Nandkumar & Pethick ( 1984). 

When differential rotation is present in a star, rotational 
energy is dissipated through frictional heating. We estimate 
the rotational energy dissipation rate as (Kippenhahn & 
Thomas [T978l Mochkovitch & Livio [T989b 

1 



rbcr^ (erg g 1 sec ^) 



(16) 



where fturb is the turbulent viscosity, which is the sum of all 
the diffusion coefficients related to the rotationally induced in- 
stabilities considered in the present study (Heger et al. l2000> . 

As discussed in Sect. 12.21 the calculation of </j^^ and (^^j 
is required for considering chemical composition effects in de- 
generate matter on the rotationally induced instabilities. We can 
make use of the thermodynamic relation: 



dlnp 
d\nP 



a + (5V , 



(17) 



where all the thermodynamic quantities on the right hand side 
are obtained from the equation of state used in our stellar evo- 
lution code (Blinnikov et al. ll996> . 

3.3. Physical assumptions 

Our CO white dwarf models start mass accretion at three differ- 
ent initial masses: 0.8, 0.9 and 1.0 Mq. Isolated white dwarfs 
are observationally found to rotate slowly (uwd Si 40 km/s ; 
Heber et al. 119971 Koester et al. 119981 Kawaler .2003i) . as also 
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Fig. 3. Evolution of central density and temperature in non- 
rotating accreting CO white dwarfs. The initial mass is 1.0 Mq 
and the accretion rate is 10"^ Mo/yr. The dotted line denotes 
the case of helium accretion, while the dashed line the case of 
carbon-oxygen accretion. The initial central temperature was 
Tc.init = 2.29 • 10* K for the helium accreting white dwarf 
and and Tcinit = 1-69 • 10* K for the CO accreting one. The 
numbers next to the filled circles denote the white dwarf mass 
in solar masses. 



predicted by stellar evolution models (Langer et al. ll999t . We 
thus consider slow rigid rotation in our initial models, and the 
surface velocity at the white dwarf equator is set to 10 km/s. 
Other physical quantities of the initial models are summarized 
in Table[T] 

As for the accretion rate, four different values are consid- 
ered: 3 • 10"^ 5 • 10"^ 7 • lO"'^ and 10"^ Mg/yr. These ac- 
cretion rates are chosen in the context of the single degenerate 
Chandrasekhar mass scenario for SNe la progenitors, in which 
steady nuclear shell burning of hydrogen or helium is assumed 
(Nomoto 1982. Li & van den Heuvel 1997. Langer et al. 2000 



Yoon & Langer l2003t . The accreted matter is carbon-oxygen 
enriched such that Xq = Xq = 0.487. In this way, we assume 
that the accreted hydrogen or helium is immediately converted 
to carbon and oxygen by shell burning. This assumption, how- 
ever, does not affect the advanced thermal evolution of the de- 
generate core of a white dwarf as demonstrated in Fig. |3l In 
this figure, the evolution of the central density and tempera- 
ture of non-rotating white dwarf which accrete matter at a rate 
of with M= 10^^ Mo/yr is shown for two cases. One case 
assumes CO accretion, the other case assumes helium accre- 
tion. In the latter case, we have followed the helium shell burn- 
ing and thereby induced thermal pulses up to AfwD= 1-2 Mq. 
Although a slight difference appears when AfwD^ 11 M©, 
which is due to the different initial temperatures as indicated in 
the figure caption, the two sequences converge when AfwD^ 
1.1 Mq due to the dominance of the compressional heating for 
the thermal structure. Therefore, we can assume safely that the 
thermal evolution of the white dwarf interior does not change 



significantly even if we neglect the shell burning, as long as 
rapid accretion is ensured. 

Since non-magnetic white dwarfs in close binary systems 
are thought to accrete matter through a Keplerian disk, the ac- 
creted matter may carry angular momentum of about the local 



Keplerian value at the white dwarf equator (cf. Durisen [19771 
Ritter I1985> . Langer et al. (JOOOj pointed out that a white 
dwarf will reach overcritical rotation well before growing to 
the Chandrasekhar limit if the white dwarf gains angular mo- 
mentum with the Keplerian value, when the assumption of rigid 
body rotation is made (see also Livio & Pringle 1998 Uenishi 
et al. 2003 1. Our preliminary calculations show that overcritical 
rotation may be reached even earlier at the white dwarf surface, 
as the angular momentum transport time scale is finite (Yoon 
& Langer 2002). Interestingly, Paczyriski J1991i) and Popham 
& Narayan t,1991 1 found a possibility for the angular momen- 
tum to be transported from the white dwarf into the disk by 
viscous effects when the white dwarf rotates at near break-up 
velocity, without preventing the continued efficient mass ac- 
cretion. Based on this picture, we limit the angular momentum 
gain such that no angular momentum is allowed to be trans- 
fered into an accreting white dwarf when its surface rotates at 
the Keplerian value at its equator: 



/ • jK if Vs < VK 

ifvs = Vk 



(18) 



where jacc denotes the specific angular momentum carried by 
the accreted matter, Vs the surface velocity of the accreting 
white dwarf at its equator, vk and jk the Keplerian value of 
the rotation velocity and the specific angular momentum at the 
white dwarf equator, respectively: vk = {GAdwD / Rwd)^^^ 
and Jk — vkRwu- The dimensionless parameter / determines 
which fraction of the Keplerian value is actually accreted in a 
white dwarf when < vk- 

As mentioned earlier, we do not consider nuclear shell 
burning in our simulations. Although this may not change the 
thermal history in the white dwarf core, the surface conditions 
can be affected by this simplification: the white dwarf envelope 
may be hotter and more extended with shell burning, which 
may lead to a stronger Eddington Sweet circulation and differ- 
ent values for jk- In order to investigate how the results are 
affected by this uncertainty, we consider three different values 
for/: 0.3, 0.5 and 1.0. 

In calibrating the chemical mixing efficiencies in massive 
stars, Heger et al. (2000i introduced the factor such that 

is replaced by /^V^, with — 0.05 giving the best re- 
production of observed enhancements of nitrogen at the stellar 
surface in the mass range 10 Mq to 20 Mq. Although this re- 
sult is deduced from massive stars and might not apply to white 
dwarfs, we keep the value of — 0.05 for most of our mod- 
els. However, in order to check the sensitivity of our results 
to this parameter, models with = 1.0 are also calculated 
(Sect. Ell. 

Table |2] gives an overview of the computed model se- 
quences. The first column gives the model sequence designa- 
tion, where the letters A, B and C denote the cases with / = 
1.0, 0.5 and 0.3 respectively, all for a fixed of 0.05. Model 
sequences with the index D have / = 1.0 and — 1.0. The 
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Table 2. Properties of the computed model sequences. The first column gives the system number. The other columns have the 
following meanings. A/jnit: initial mass, M: accretion rate, /: fraction of jk of the accreted matter (see Ea.ll8>. : efficiency 
parameter of the ji gradient, Erot'- rotational energy dissipation by friction, A/q.io- white dwarf mass when the ratio of the 
rotational energy to the gravitational potential energy (T/W) reaches 0.10. The rest physical quantities are also estimated at this 
point, i.e., when T/W — 0.10. Pc,o.io and Tc^o.io '■ central density and temperature. i?o.io: radius of white dwarf defined on the 
sphere of the volume of the equipotential surface. Jo.io ■ total angular momentum, ilo.io: moment of inertia weighted mean of 
angular velocity. 



No. 


Minit 


M 


f 


U 


Eiot 


Mo.w 


Pc.O.lO 


7c, 0.10 


.Ro.io 


Jo. 10 


f^o.io 






10-' 










10'' 


10" 


0.01 


10^" 






Mq 


Me/yr 








Mq 


g/cm^ 


K 


R© 


erg s 


rad/s 


Al 


0.8 


3.0 


1.0 


0.05 


Yes 


1.16 


0.29 


0.92 


1.22 


1.10 


0.68 


A2 


0.8 


5.0 


1.0 


0.05 


Yes 


1.18 


0.32 


0.97 


1.26 


1.13 


0.70 


A3 


0.8 


7.0 


1.0 


0.05 


Yes 


1.19 


0.33 


1.01 


1.34 


1.14 


0.69 


A4 


0.8 


10.0 


1.0 


0.05 


Yes 


1.21 


0.35 


1.07 


1.39 


1.17 


0.71 


A5 


0.9 


3.0 


1.0 


0.05 


Yes 


1.28 


0.58 


1.02 


1.02 


1.18 


0.95 


A6 


0.9 


5.0 


1.0 


0.05 


Yes 


1.30 


0.64 


1.11 


1.03 


1.21 


0.99 


A7 


0.9 


7.0 


1.0 


0.05 


Yes 


1.32 


0.69 


1.17 


1.03 


1.23 


1.02 


A8 


0.9 


10.0 


1.0 


0.05 


Yes 


1.34 


0.75 


1.23 


1.05 


1.25 


1.06 


A9 


1.0 


3.0 


1.0 


0.05 


Yes 


1.41 


1.30 


1.14 


0.82 


1.26 


1.40 


AlO 


1.0 


5.0 


1.0 


0.05 


Yes 


1.42 


1.39 


1.22 


0.83 


1.25 


1.43 


All 


1.0 


7.0 


1.0 


0.05 


Yes 


1.43 


1.52 


1.28 


0.82 


1.27 


1.49 


A12 


1.0 


10.0 


1.0 


0.05 


Yes 


1.45 


1.67 


1.36 


0.83 


1.28 


1.56 


Bl 


0.8 


3.0 


0.5 


0.05 


Yes 


1.21 


0.36 


0.96 


1.17 


1.17 


0.76 


B2 


0.8 


5.0 


0.5 


0.05 


Yes 


1.25 


0.42 


1.04 


1.20 


1.22 


0.81 


B3 


0.8 


7.0 


0.5 


0.05 


Yes 


1.26 


0.44 


1.09 


1.24 


1.23 


0.82 


B4 


0.8 


10.0 


0.5 


0.05 


Yes 


1.30 


0.52 


1.16 


1.23 


1.27 


0.89 


B5 


0.9 


3.0 


0.5 


0.05 


Yes 


1.32 


0.69 


1.03 


0.97 


1.24 


1.04 


B6 


0.9 


5.0 


0.5 


0.05 


Yes 


1.35 


0.79 


1.12 


0.97 


1.27 


1.11 


B7 


0.9 


7.0 


0.5 


0.05 


Yes 


1.38 


0.90 


1.12 


0.96 


1.29 


1.18 


B8 


0.9 


10.0 


0.5 


0.05 


Yes 


1.41 


1.04 


1.27 


0.97 


1.32 


1.26 


B9 


1.0 


3.0 


0.5 


0.05 


Yes 


1.45 


1.61 


1.18 


0.78 


1.30 


1.56 


BIO 


1.0 


5.0 


0.5 


0.05 


Yes 


1.46 


1.76 


1.27 


0.78 


1.30 


1.62 


Bll 


1.0 


7.0 


0.5 


0.05 


Yes 


1.48 


2.00 


1.34 


0.78 


1.31 


1.72 


B12 


1.0 


10.0 


0.5 


0.05 


Yes 


1.51 


2.33 


1.44 


0.77 


1.32 


1.84 


CI 


0.5 


J.U 


U.3 


U.Uj 


Yes 


1.D3 


1.03 


1 1 "7 
1.1 / 


U. / 1 


1.45 


1.3 / 


C2 


0.8 


5.0 


0.3 


0.05 


Yes 


1.50 


1.25 


1.19 


0.84 


1.46 


1.41 


C3 


0.8 


7.0 


0.3 


0.05 


Yes 


1.46 


1.05 


1.20 


0.90 


1.44 


1.29 


C4 


0.8 


10.0 


0.3 


0.05 


Yes 


1.42 


0.85 


1.21 


1.01 


1.40 


1.15 


C9 


1.0 


3.0 


0.3 


0.05 


Yes 


1.73 


19.0 


2.16 


0.42 


1.25 


5.03 


CIO 


1.0 


5.0 


0.3 


0.05 


Yes 


1.72 


16.2 


2.04 


0.44 


1.27 


4.65 


CU 


1.0 


7.0 


0.3 


0.05 


Yes 


1.72 


14.3 


2.06 


0.46 


1.29 


4.38 


C12 


1.0 


7.0 


0.3 


0.05 


Yes 


1.71 


12.1 


2.01 


0.49 


1.31 


4.04 


D2 


0.8 


5.0 


1.0 


1.00 


Yes 


1.18 


0.32 


0.97 


1.27 


1.13 


0.69 


D6 


0.9 


5.0 


1.0 


1.00 


Yes 


1.30 


0.64 


1.11 


1.03 


1.21 


0.99 


T2 


0.8 


5.0 


1.0 


0.05 


No 


1.12 


0.27 


0.95 


1.24 


1.02 


0.63 


T6 


0.9 


5.0 


1.0 


0.05 


No 


1.26 


0.56 


1.03 


0.97 


1.15 


0.93 


TIO 


1.0 


5.0 


1.0 


0.05 


No 


1.40 


1.24 


1.14 


0.81 


1.24 


1.37 



index T indicates test sequences where the effects of rotational 
energy dissipation (Eq.O are neglected. 

4. Results 

In this section, we present the results of our simulations, which 
are summarized in Tables |2] and |3] In these tables, key prop- 
erties of the white dwarf models are given for three different 
epochs, i.e., when the ratio of the rotational energy to the gravi- 
tational potential energy Emt /\W\ reaches 0. 10, 0. 14 and 0.18. 



The implications of these numbers are discussed in Sect. [S] In 
the following subsections, after discussing the thermal evolu- 
tion of our models, we focus on the processes of the angular 
momentum transport in the interior of accreting white dwarfs 
and its consequences for the white dwarf structure. 

4. 1. Thermal evolution 

The thermal evolution of non-rotating accreting white 
dwarfs has been studied by many authors (e.g. Iben 119821 
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Table 3. Continued from Table|2] The symbols in the columns have the same meaning with those in Table|3l but when T/W — 
0.14 and T/W = 0.18 for the indices 0.14 and 0.18, respectively. 



No. 


Mn 1 1 


pc,0. 14 


c,0. 14 


— P 

14 


— J 

'Jo. 14 


~o 1 

i '0. 14 


Mn iH 


pc.0. 18 


J- c,0. 18 


— P 

^0.18 


— J 

"0.18 


" ^0.18 






10** 


10** 


0.01 


10^" 






10*^ 


10** 


0.01 


10^" 








&/ '-"i 


K 




pro 


r3.d/s 


iV10 


&/ "-"1 


K 




pre 


rsd/s 


Al 


1 .J 1 


0.42 


04 


1 10 


1 S4 


Q7 


1 4'i 


'Sfi 


06 


1 01 

1 .U i 


1 Q7 


i .z,u 


A2 


1 "^4. 
1 .j^ 


n 47 


1 01 


1.11 


1 


1 07 


1.47 




1 04 


1 01 

1 .U i 






A3 


1.35 


0.49 


1.06 


1.13 


1.59 


1.03 


1.48 


0.66 


1.09 


1.03 


2.00 


1.34 


A4 




n S4 


1.12 


1 1 


1.63 


1 07 


1 'SO 


7^ 


1 1 s 


1 04 


9 0^ 


1 40 


A5 


1.42 


0.78 


1.03 


0.93 


1.62 


1.30 


1.55 


1.04 


1.74 


0.87 


2.02 


1.68 


A6 


1.44 


0.87 


1.11 


0.93 


1.63 


1.36 


1.56 


1.15 


1.15 


0.87 


2.01 


1.74 


A7 


1.46 


0.97 


1.17 


0.93 


1.65 


1.42 


1.58 


1.28 


1.21 


0.85 


2.03 


1.83 


A8 


1.48 


1.07 


1.23 


0.93 


1.67 


1.48 


1.60 


1.42 


1.28 


0.85 


2.04 


1.92 


A9 


1.55 


1.78 


1.18 


0.76 


1.66 


1.90 


1.69 


2.56 


1.28 


0.69 


2.03 


2.54 


AlO 


1.55 


1.88 


1.26 


0.77 


1.65 


1.94 


1.67 


2.55 


1.33 


0.71 


2.01 


2.52 


All 


1.57 


2.08 


1.33 


0.76 


1.66 


2.03 


1.69 


2.82 


1.40 


0.70 


2.00 


2.63 


A12 


1.58 


2.31 


1.40 


0.75 


1.66 


2.12 


1.70 


3.13 


1.49 


0.70 


2.00 


2.75 


Bl 


1 "^fi 
1 . 


0.51 


U.7U 


1 Oft 


1 (SI 


1 07 


1 'SO 


71 

\J. 1 L 


1 00 
1 .uu 


07 
\j.y 1 


9 0^ 


1.41 


B2 


1 4.0 


u.oz. 


1 ns 

1 .UJ 


1 OS 


1 fif< 


1 1(S 

1 . lU 


1 'S4 





1 00 
1 .\jy 


06 


2.07 


1 

L.J J 


ijj 


1.41 




1 10 


1 05 


1 (S7 




1 'S'S 

L.J J 


so 


1.14 


06 


2.07 


1 SS 

L.J J 


B4 




78 
u. / o 


1.17 


1 06 


1.71 


1.28 


1 'S8 

L.JO 


1 07 
1 .u / 


1.22 


OS 


9 00 


1 60 


B5 


1.46 


0.95 


1.06 


0.90 


1.66 


1.43 


1.60 


1.32 


1.12 


0.82 


2.06 


1.88 


B6 


1.49 


1.12 


1.15 


0.89 


1.68 


1.54 


1.63 


1.56 


1.22 


0.81 


2.07 


2.02 


B7 


1.52 


1.28 


1.22 


0.86 


1.70 


1.63 


1.65 


1.79 


1.29 


0.81 


2.07 


2.15 


B8 


1.55 


1.52 


1.31 


0.87 


1.72 


1.76 


1.67 


2.12 


1.39 


0.77 


2.07 


2.32 


B9 


1.60 


2.38 


1.36 


0.71 


1.70 


2.19 


1.75 


3.94 


1.40 


0.62 


2.05 


3.13 


BIO 


1.60 


2.52 


1.34 


0.72 


1.68 


2.23 


1.74 


3.81 


1.46 


0.64 


2.02 


3.05 


Bll 


1.62 


2.90 


1.42 


0.70 


1.68 


2.38 


1.75 


4.32 


1.54 


0.64 


2.01 


3.22 


B12 


1.64 


3.43 


1.53 


0.68 


1.68 


2.56 


1.75 


4.91 


1.65 


0.61 


1.98 


3.42 


CI 


1.74 


3.36 


1.36 


0.64 


1.87 


2.62 


1.87 


6.30 


1.57 


0.54 


2.13 


3.99 


C2 


1.71 


2.66 


1.38 


0.69 


1.89 


2.35 


1.86 


5.02 


1.58 


0.59 


2.17 


3.58 


C3 


1.68 


2.18 


1.37 


0.74 


1.89 


2.13 


1.84 


4.13 


1.57 


0.63 


2.20 


3.27 


C4 


1.64 


1.69 


1.35 


0.81 


1.88 


1.88 


1.81 


3.20 


1.54 


0.69 


2.23 


2.89 


D2 


1.34 


0.47 


1.01 


1.11 


1.58 


1.02 


1.47 


0.63 


1.04 


1.01 


1.99 


1.32 


D6 


1.44 


0.87 


1.11 


0.94 


1.63 


1.36 


1.56 


1.15 


1.15 


0.87 


2.01 


1.74 


T2 


1.26 


0.37 


0.95 


1.12 


1.45 


0.90 


1.38 


0.47 


0.95 


1.05 


1.85 


1.14 


T6 


1.39 


0.72 


1.02 


0.94 


1.57 


1.25 


1.51 


0.72 


1.02 


0.94 


1.57 


1.25 


TIO 


1.52 


1.57 


1.15 


0.76 


1.63 


1.79 


1.64 


2.06 


1.21 


0.71 


2.00 


2.29 



Nomoto [T9821 Sion [T9951 Townsley & Bildsten l2002l . 
Piersanti et al. ( 120031 compared rigidly rotating white dwarf 
models with non-rotating ones for CO-accretion at various ac- 
cretion rates. They found that a rotating white dwarf is gener- 
ally cooler than its non-rotating counterpart, which is a natural 
consequence of the lifting effect of the centrifugal force. This 
conclusion holds also in our differentially rotating white dwarf 
models (Fig.|4}. However, as we shall see in Sect. 14. 21 the ther- 
mal structure in an accreting white dwarf is found to have an in- 
teresting consequence for the redistribution of angular momen- 
tum, as it affects the stability criterion for the dynamical shear 
instability. Therefore, before discussing the angular momentum 
transport, we examine the evolution of the thermal structure in 
our white dwarf models. 

Fig. |4] displays the temperature profiles of selected white 
dwarf models of different sequences. The non-rotating mod- 
els (Fig. 1^) show a rapid increase of the central tempera- 
ture, from 10^ K when Mwd = 0.96 Mq to 2.7 • 10** K 
when Mwd = 1-37 M©. In the corresponding rotating mod- 



els (Fig. 13) & c), the change in the central temperature is not 
significant in the considered period. 

A comparison of Fig. and Fig. |3 shows the effects 
of rotational energy dissipation (-Erot) described by Eq. 
In model sequence T6 (Fig. where i^iot is neglected, the 
white dwarf models are cooler than in the corresponding mod- 
els with E-cot considered (sequence A6, Fig.|3). As shown in 
Fig. Is] the rotational energy dissipation rate can be as high as 
10'^ . . . lO'' erg/g/s in the degenerate core, and even higher 
close to the surface. The integrated energy dissipation rate 
Erot — J ^dissdM^ reaches several times 10^ L0 in these mod- 
els. This results in a heating of the region with a strong shear 
motion (cf. Fig.|6lbelow), with a consequence that the temper- 
ature maxima when A/wd = 0.96, 1.22 and 1.37 Mq in models 
of sequence A6 are higher than in the models of sequence T6. 
The accretion induced heating is a sensitive function of the ac- 
cretion rate: comparison of Fig.|3; with Fig.^ indicates that 
the white dwarf becomes significantly hotter with a higher ac- 
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Fig. 4. Temperature as function of the mass coordinate in white dwarf models when Mwd = 0.96, 1.22 and 1.37 M©. (a) non- 
rotating case with Minit = 0.9 Mq and M = 5 • lO^*^, (b) model sequence T6, (c) model sequence A6, (d) model sequence 
A6. 



cretion rate. The effects of the accretion rate and of £^iot for the 
angular momentum transport are discussed in Sect. 14.3.31 

Note that in all white dwarf models, there exists an abso- 
lute temperature maximum in the outer layers of the degenerate 
core. Below this temperature peak, the temperature gradient V 
becomes negative and large. This produces a strong buoyancy 
force in this region, since the Brunt- Vaisala frequency is pro- 
portional to Vad — V. The temperature peak moves outward as 
the white dwarf accretes more mass. This leads to changes in 
the local thermodynamic condition which plays a key role for 
the angular momentum transport from the outer envelope into 
the inner core, as shown in the next section. 

4.2. White dwarf spin and angular momentum 
transport 

Our white dwarf models are spun up as they gain angular mo- 
mentum from the accreted matter, as described in Sect. 13.31 
Fig- El shows the evolution of the angular velocity profiles in 
model sequences A2, A6, AlO and D6. Most strikingly, all 
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Fig. 5. Rotational energy dissipation rate (see Eq. I16> as func- 
tion of the mass coordinate, when Mwd = 0.98, 1.22 and 1.37 
M0 in model sequence A6. 
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Fig. 6. Angular velocity as function of the mass coordinate at different white dwarf masses. Panels (a), (b), (c) & (d) give results 
for model sequences A2, A6, AlO and D6, respectively (see Table|2j. 



2.0 



1.5 



B 1.0 



0.0 







, 




_ Mtait = 0.9 Mo 


1 




" dM/dl = 5-10-'' Mo/yr 










a : MwD = 1.26 M© 










b : MwD = 1.37 Me 


a / 










^ b 








^1 , , 1 , , , 1 , , , 1 , , 


1 < , < 1 < < 




il 1 iT 



0.0 0.2 0.4 0.6 O.B 1.0 1.2 1.4 
Mr/Mo 



Fig. 7. Shear factor a (= doj/dhir) as function of the 
mass coordinate in two different white dwarf models at 
-^WD = 1-26 M0 and 1.37 M0 from model sequence A6. The 
dashed lines denote the threshold value of a for the dynamical 
shear instability (i.e., CTDSLcrit — (A^Vi?i,c)^/^ see Eq.EJ. 



white dwarf models rotate differentially, as predicted from the 
discussions in Sect. |2l Note also that in every white dwarf 
model in Fig.|6l a maximum angular velocity occurs, such that 
(7 > in the inner core and a < in the outer layers. This max- 
imum comes into existence soon after the onset of mass accre- 
tion because the slowly rotating inner part contracts faster than 
the rapidly rotating surface layers as the total mass increases. 

Since angular momentum is transported in the direction of 
decreasing angular velocity, this peak in lu serves as a bottle 
neck for the angular momentum transport from the outer en- 
velope into the core. It prevents the outer envelope from slow- 
ing down efficiently by inward angular momentum transport. 
As a result, the surface remains to rotate close to the critical 
value throughout the evolution, severely limiting the angular 
momentum gain from the accreted matter. I.e., about 60 % of 
the angular momentum of the accreted matter is rejected by the 
condition posed in Eq. ^] and only about 40 % is actually re- 
tained by the time the white dwarf mass reaches 1.4 Mq, in all 
model sequences (see Fig.ll llbelow). 

Fig.0shows the shear factor a as function of the mass co- 
ordinate at two different evolutionary epochs (AfwD = 1-26 & 
1.37 Mq) in model sequence A6. The dashed line gives the 
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threshold value of the shear factor (fiDSi.crit)- Note that the 
two, i.e., a and (TDSi,crit converge remarkably well in the white 
dwarf interior, i.e., in M-^ < 1.1 Mq when Af-yvD = 1-26 Mq 
and Mt- < 1.2 Mq when Af^vD = 1-37 Mq. This confirms 
the conclusion given in Sect. |2] that the degenerate core of an 
accreting white dwarf will rotate differentially with the shear 
strength near the threshold value for the dynamical instability. 
As already discussed in Sect. 12.21 this is because any strong 
shear motion with a > CDSi.crit can not be retained for a long 
time and should decay to CTDSi.crit quickly via the dynamical 
shear instability, and because further angular momentum trans- 
port by other mechanisms requires a longer time scale com- 
pared to the accretion time scale. 

Fig. |8l depicts the situation in more detail. The right hand 
panel of this figure shows the various diffusion coefficients as 
function of the mass coordinate in the white dwarf models of 
sequence A6 when M^d = 1-26 and 1.37 Mq. This figure 
shows that the white dwarf consists of an inner dynamical shear 
unstable region and an outer region dominated by the secular 
shear instability and Eddington Sweet circulations. Let us de- 
fine A/j-.DSi as the mass coordinate at the point which divides 
the white dwarf into these two regions. A/i,dsi changes with 
time not only due to the change of the shear strength, but also 
due to the evolution of the thermodynamic properties as ex- 
plained below. 

As shown in Sect. 14. II the accretion induced heating results 
in a temperature maximum in the outer region of an accreting 
white dwarf. A steep negative temperature gradient (i.e., V < 
0) thus appears just below this temperature peak. The buoyancy 
force is enhanced in this region, which leads to stability against 
the dynamical shear instability (See Eqs.|3}|6}. Therefore, the 
location of Af^ dsi is below the region containing the strong 
temperature gradient, as demonstrated in Fig.|8] 

According to Fig. |8] A/i dsi moves outward as the white 
dwarf mass increases (See also Fig. and Fig. |9jl. This out- 
ward shift of A/i,DSi can be understood as follows. As shown 
in Fig. the threshold value crDSi,ciit becomes smaller near 
Afr,DSi as the white dwarf mass increases. Two processes con- 
tribute to this effect. First, the degeneracy at a given Air be- 
comes stronger as the white dwarf mass increases, reducing 
the buoyancy force. Second, the region with strong tempera- 
ture stratification with V < moves outward as indicated in 
Fig Et and Fig.lSj;. Consequently the dynamical shear unstable 
region gets extended outward with time as shown in Fig.|8land 
also in Fig.|9l The outward shift of M^^usi allows the angular 
momentum from the outer layers to be transported inward, and 
thus leads to the outward shift of the position of the maximum 
angular velocity as shown in Fig.|9] 

Fig. ^1 shows how much angular momentum is actually 
transported into the white dwarf interior in model sequence 
A6. Here J{M-c) — J^'' j{m)dm is the integrated angular 
momentum. In a rigidly rotating homogeneous body, we have 
J[M-c)/Mr^'^ = const., and therefore any deviation from the 
zero gradient in J{M,-)/Mr^'^ may serve as a measure of the 
degree of differential rotation. A comparison of different evo- 
lutionary epochs demonstrates the flow of angular momentum 
through the mass shells, as J(A/i.) and J{M-c)/M^^^ would re- 
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Fig. 9. Angular velocity profiles given as function of the mass 
coordinate at 17 different evolutionary epochs of sequence A6: 
1.26, 1.28, 1.30, 1.31, 1.33, 1.34, 1.36, 1.37, 1.38, 1.40, 1.42, 
1.44, 1.45, 1.47, 1.48, 1.50, 1.52 Mq, from the bottom to the 
top. The filled circles designate the mass coordinate A/rm,DSi> 
below which the white dwarf interior is dynamical shear unsta- 
ble. 

main constant at a given mass shell throughout the evolution 
if no angular momentum were transported through this shell. 
Note that all the white dwarf models in this figure have a sim- 
ilar amount of angular momentum (J ~ 10''^" erg s) at around 
A/f ~ 1.2 Mq, which is close to the angular velocity peaks 
shown in Fig. |9] The angular momentum transport from this 
point into the further interior is shown to be efficient, due to 
the dynamical shear instability as described above. 

We summarize the main conclusions from this section as 
follows, a) We find accreting white dwarf models to rotate dif- 
ferentially throughout their evolution, b) The dynamical shear 
instability is the most important mechanism for angular mo- 
mentum transport in the highly degenerate core of our white 
dwarf models, c) The angular momentum gain from the ac- 
creted matter and the spin-up of our accreting white dwarfs is 
closely related to their restructuring and thermal evolution as 
their mass increases. 

4.3. Influence of physical assumptions 

Having understood the detailed processes of the angular mo- 
mentum transport in an accreting white dwarf, let us consider 
the influence of the different initial conditions and physical as- 
sumptions on the evolution of accreting white dwarfs. 

4.3.1 . Influence of the initial white dwarf mass 

Tables |2] and |3] give the white dwarf mass when the ratio 
of the rotational energy to the gravitational potential energy 
(i?iot/| W^l) in each model sequence reaches 0.1, 0.14 and 0.18 
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Fig. 8. (a) Thermodynamic quantities as function of the mass coordinate, in the white dwarf model of sequence A6 when 
A/wD = 1-26 Vad (solid line), V (dot-dashed line), /p(<PcVp^ + (pi\/^^)/6 where = 0.05 (dashed line) and temperature 
(dotted line). The region where the white dwarf is unstable to the dynamical shear instability is hatched, (b) Diffusion coefficients 
for rotationally induced hydrodynamic instabilities, as function of the mass coordinate, in the white dwarf model of sequence A6 
at AfwD = 1 -26 Mq . The solid line denotes the diffusion coefficient for the dynamical shear instability, the dashed line, the dotted 
Une and the dashed-dotted lines are for the secular shear instability, the Eddington Sweet circulation, and the GSF instability 
respectively, (c) Same as in (a) but when Mwd = 1-27 Mq (d) Same as in (b) but when AfwD = 1-37 Mq 



(Mo,i, Afo.i4 and Mq.is)- It is believed that a rapidly rotating 
white dwarf becomes secularly unstable when Erot / 1 VF | ~ 0. 14 
(e.g. Durisen & Imamura'1981V Imamura et al (^1995) found 
that this critical value can be lowered to about 0. 1 for strong 
differential rotation (see detailed discussion in Sect.|5|l. 



The tables show that the lower the initial mass, the higher 
the values of Mq iq and M0.14 become. For example, A/0.14 
increases from 1 .34 Mq to 1 .55 Mq with a change of the initial 
mass from 0.8 Mq to 1.0 Mq in the case of / = 1 (sequences 
A2 and A 10). The total angular momentum also increases at 
a given i^rot/lW^I with higher initial mass. This tendency is 
simply due to the fact that more angular momentum gain is 
necessary to reach a certain amount of E'rot/lW^I if a white 
dwarf is more massive. 



4.3.2. Influence of the angular momentum gain 
parameter / 

As indicated in Table. |2land|3] model sequences with / = 1.0 
and those with / = 0.5 show some differences in the results. 
The adoption of / = 0.5 yields a larger white dwarf mass at a 
given i?rot/|W^|, by about 0.05 - 0.06 Mq, which is a natural 
consequence of less angular momentum being accreted per unit 
time (Eg. list . This effect becomes more prominent with / = 
0.3, and Afo.i increases by about 24 % and 17 % for A/jnit = 
0.8 Mq and Afinit = 1.0 Mq respectively, compared to the case 
of / — 1. Furthermore, the models with A/init =1.0 Mq and 
/ = 0.3 reach central carbon ignition when £^rot/|M^| reaches 
about 0.11. 

Interestingly, the helium accreting white dwarf models by 
Yoon et al. (.2004J . where / = 1 is adopted, show similar 
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Fig. 10. Integrated angular momentum J(Afi) = 

J^' j{m)dm divided by Mr^'^, as function of the mass 
coordinate at five different evolutionary epochs of sequence 
A6: MwD = 1-21, 1.30, 1.37, 1.44 and 1.52 Mq. The thin 
contour lines denote levels of constant J in logarithmic scale, 
labeled with log( J/10^"ergs). 

i'rot/|W^| values at a given mass as the present CO accreting 
white dwarf models with / = 0.3. As mentioned in Sect. 13.21 
this difference is mainly due to the fact that the angular mo- 
mentum transport efficiency in the outermost layers is affected 
by the energy generation in helium burning shell. In particular, 
the Eddington Sweet circulation becomes more efficient due to 
shell burning, resulting in a more efficient outward angular mo- 
mentum transport in the non-degenerate envelope (cf. Fig.|6}. 
This causes a more severe restriction of the angular momen- 
tum gain from the accreted matter, due to the condition posed 
by Eq.^] This implies that the history of the angular momen- 
tum gain may also be different for hydrogen accreting cases. 

However, these ambiguities concern only the actual amount 
of the angular momentum gain from the accreted matter, but do 
not affect the history of the angular momentum redistribution in 
the degenerate core, where differential rotation persists during 
the mass accretion phase (Sect. IZ21 and I4.2l l. In fact, despite 
big differences in / all model sequences show the same re- 
markable feature the carbon ignition is not reached even when 
M^Y^ > 1.4 Mq, due to differential rotation. 

4.3.3. Influence of accretion rate and dissipation 

Tables 121 and |3l show that with a given initial mass, the white 
dwarf mass at at given i^rot/lM^I increases for higher accre- 
tion rate. Fig. ^3 gives the accumulated angular momentum 
(A Jaccictcd) in the white dwarf, as well as the accumulated 
rejected angular momentum (A Jrojcctod) according to Eg. 1181 
as function of the white dwarf mass, for sequences with A/jnit 
= 0.9 Mq, and for three different accretion rates as indicated 
in the figure. It is shown that the higher the accretion rate, the 
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Fig. 11. The evolution of the accumulated angular momentum 
in model sequences A5, A6 and AlO, given as function of the 
white dwarf mass. The thin lines denote the accumulated total 
angular momentum in the white dwarf models. The thick lines 
give the accumulated rejected angular momentum by Eg. 1181 

more angular momentum is rejected and the less angular mo- 
mentum is accreted, which is the reason for the higher white 
dwarf mass at a given i?rot/|W^| for a higher accretion rate. 

This accretion rate dependence can be explained by two 
factors. First, with a lower accretion rate, a white dwarf gains 
a smaller amount of angular momentum per time and thus 
has more time to transport angular momentum into the white 
dwarf interior Second, as already pointed out in Sect. 14.11 a 
higher accretion rate results in a stronger accretion induced 
heating, which leads to higher temperatures inside the white 
dwarf (Fig.|4}. This reduces the degeneracy in the white dwarf 
interior, and the buoyancy force becomes accordingly stronger. 
This thermal effect changes the stability condition for the dy- 
namical shear instability in the white dwarf interior: the higher 
the accretion rate, the less susceptible to the dynamical shear 
instability it is, which in turn limits the angular momentum 
transfer from the outer envelope into the interior more severely. 

We can also understand why the white dwarf mass at a 
given E'rot/lW^I is smaller in case E'l-ot is neglected (Table |2] 
and O in terms of different thermal structures: the white 
dwarf temperature becomes significantly lower without E^ot 
as shown in Fig.|3 resulting in a lower buoyancy force in the 
white dwarf interior and thus a more efficient angular momen- 
tum transport. 

4.3.4. Influence of chemical gradients 

In the model sequences discussed in the previous sections, the 
effect of (U-gradients is significantly reduced by using = 
0.05, as described in Sect. 13.31 In order to understand the 
importance of chemical gradients for the angular momentum 
evolution, two model sequences (D2 and D6) are computed 
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with = 1.0. The results are presented in Table |2] and |3] 
Interestingly, it is found that the adoption of different has 
hardly any effect. For instance, we have exactly the same val- 
ues for Afo.i, Afo.i4 and A/q.is in the two corresponding model 
sequences A2 and D2, as well as in A6 and D6. 

Fig. I12h shows the chemical structure in two white dwarf 
models of sequence D6, in the initial model (0.9 Mq) and 
when A/wD = 1-37 Mq. The model at 1.37 Mq shows that al- 
though the rotationally induced chemical mixing smoothes out 
the chemical structure significantly, a strong stratification in the 
chemical composition persists in the range Afr = 0.4— 0.8 Mq. 
Nevertheless, its contribution to the total buoyancy force turns 
out to be too small to suppress the dynamical shear instability, 
as implied in Fig. 112b. The term ((^oV^^ + ipi\/^^)/6 has a 
maximum value of about 0.1, which is only 25 % of Vad- 

In conclusion, the effect of chemical gradients on the angu- 
lar momentum transport is negligible in the white dwarf inte- 
rior, unlike in non-degenerate stars. 

4.3.5. Limitations of our 1 D approximation 

As pointed out in Sect. 13. 31 our one-dimensional description of 
rotation works accurately up to about 60% of Keplerian rota- 
tion. In our white dwarf models, however, the fast rotating outer 
layers exceed this limit (cf. Fig.ll3t. For instance, in model se- 
quence A6, the outer 19 % in mass rotate faster than 60 % of 
the Keplerian value when i?rot/|M^| =0.1. This fast rotating 
region increases to 27 % in mass when i?rot/|W^| reaches 0.14. 
This means that our numerical models underestimate the effect 
of the centrifugal force on the white dwarf structure. 

The accreting white dwarf models by Durisen fl977^ may 
be most appropriate to evaluate the uncertainty due to the un- 
derestimation of the centrifugal force in the outer envelope. 
In his two dimensional models, the inner dense core with 
p > O.lpc remains close to spherical, and the outer layers with 
p < 0.1 pc start deviating from the spherical symmetry, becom- 
ing toroidal. Consequently, the outer envelope is considerably 
extended, making the ratio of the polar radius to the equatorial 
one (Rp/Rc) as small as 0.4 when £'rot/|W^| reaches 0.1. 

In our models, the inner slowly rotating core (uj < 
O.GwKcpior) has density p > O.lpc in general, and therefore 
we may conclude that the inner core is accurately described 
in our calculations, although it could be implicitly affected by 
the underestimation of the centrifugal force in the outer enve- 
lope. However, the major qualitative conclusions of our study 
are not affected by this uncertainty: accreting white dwarfs 
will rotate differentially and can grow beyond the canonical 
Chandrasekhar mass without suffering central carbon ignition, 
unless they lose angular momentum through secular instabili- 
ties (Sect.|5}. 

5. On the final fate of rapidly rotating massive 
white dwarfs 

A remarkable feature in our models is that the white dwarfs 
can not reach central carbon ignition even when their mass be- 
comes significantly larger than the canonical Chandrasekhar 




Fig. 12. (a) Mass fraction of carbon as function of the mass co- 
ordinate in two white dwarf models of sequence D6, at AfyvD 
= 0.9 Mq (initial model, dashed line) and 1.37 Mq (solid line), 
(b) Thermodynamic quantities as function of the mass coor- 
dinate for the same models: Vad (solid line), V (dot-dashed 
line), (iy9oV^^ + iy9iVpj)/(5 (dashed line) and temperature (dot- 
ted line). The region where the white dwarf is unstable to the 
dynamical shear instability is hatched. 

mass of 1.4 Mq (Table |2] and |3}- For example, in model se- 
quence A12, the central density at A/wd = 1-70 Mq is only 
3.13 X 10* g/cm"^ (see Table|3, which is still far from the car- 
bon ignition density (a few 10^ g/cm'^ at T > 10* K). This 
is in good agreement with the conclusion of previous studies 
of rotating white dwarfs (cf. Sect. 13. 1> that a white dwarf can 
be dynamically stable for masses up to ~ 4.0 Mq if it rotates 
differentially. 

If a rapidly rotating massive white dwarf would not lose 
angular momentum, it could not produce a Type la super- 
nova. However, it is well known that rapidly rotating com- 
pact stars becomes secularly unstable to non-axisymmetric 
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Fig. 13. Angular velocity normalized to the local Keplerian 
value, as function of the mass coordinate in white dwarf models 
of sequence A6, when E,ot/\W\ = 0.04, 0.07, 0.10 and 0.14 



perturbations due to the gravitational wave radiation reaction 
(Chandrasekhar [T970l Friedman & Schutz[T978j, which is of- 
ten named CFS (Chandrasekhar, Friedman and Schutz) insta- 
bility. Two representative modes are believed to be the most 
important ones for this instability. One is the bar-mode, i.e., 
the /-mode with ^ = m = 2, of which the restoring force 
is the buoyancy force. Here, I and ni denote the nodal num- 
bers of the spherical harmonics. Although Friedman & Schutz 
J 197 8') found that /-modes with a higher m become more sus- 
ceptible to the CFS instability compared to the case m — 2, 
their growth time becomes usually too long to be of astrophys- 
ical interest (e.g. Shapiro & Teukolskv ll983> . Hereafter, we 
will refer the CFS instability with m = 2 as the bar-mode 
instability. Recently it has been found that the CFS instabil- 
ity can also excite the 7--mode, of which the restoring force is 
the Coriolis force (hereafter, r-mode instability, see Andersson 
& Kokkotas 2001 for a review). Here again, the mode with 
m = ^ = 2 is most relevant because it gives the smallest time 
scale for for the growth of the instability. 

In the following, we discuss the importance of these insta- 
bilities for our white dwarf models, and derive implications for 
the final fate of rapidly rotating massive white dwarfs. 

5.1. Bar-mode instability 

The bar- mode instability can operate when Srot/IM^I exceeds 
a certain critical value [(S'rot/lW^Dbar, crit] in a rotating star, as 
studied by many authors (Ostriker & Tassoul il969. Ostriker & 
Bodenheimer [TW3l Durisen[T975b|[T977l Bardeen et al. [T977l 
Durisen & Imamura ' 1981t . Although (i?rot/|W^|)bar, crit is 
found to be about 0.14 for a wide range of rotation laws 
and equations of states (e.g. Durisen I1975bl Karino & 
Eriguchi '20021, Imamura et al. ( 19951 showed that it tends to 
decrease for strong differential rotation. In particular, we note 
that the rotation law in our models bears a similarity with one 
of their rotating polytrope models with n' = oo (see Imamura 
et al. for the definition of n'), in which the spin rate shows a 



maximum where the gradient in u changes its sign, as in our 
models (Fig.|6}. In these models, (i?rot/|W^|)bar, crit turns out 
to be as small as 0.09, which is significantly smaller than the 
canonical value of 0. 14. 

According to Chandrasekhar J1970t and Friedman & 
Schutz 119751 . the growth time of the bar mode instability in 
Maclaurin spheroids is given by 
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(19) 

for < E^.ot/\W\ - {E^-ot/\W\)c < 1, where fl is the mean 
angular velocity and R is the radius. 

Since we have Erot ^ k^MR^n^ and \W\ ~ GAP/R 
(where k is the dimensionless radius of gyration), the growth 
time of the instability can be rewritten as (cf. Friedman and 
Schutz [T9751 Hayashi et al. [T998l 
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where = 0.4 has been used. Although our white dwarf mod- 
els deviate significantly from the Maclaurin spheroids, this ap- 
proximation may give an order of magnitude estimate for the 
growth time scale of the bar-mode instability. 

Fig.ll4lshows Tbar as function of Srot/I for four differ- 
ent masses as indicated in the figure caption. Here, R = 0.01 
R0 is assumed since the mean radii of our white dwarf mod- 
els do not differ much from 0.01 R© for all i?rot/|W^| values. 
Two different values for (i^rot/lW^Dbar. crit, i.e., 0.1 and 0.14 
are considered in the figure. The growth time scale of the bar- 
mode instability is a sensitive function of E'rot/lW^I, while it 
does not change much for different white dwarf masses: In 
case of (£^rot/|M^|)bar, crit = 0.10, Tbar is as large as 10^° yr 
at Erotl\W\c:i (-Erot/liy Dbar, crit, but it drops to 10^ yr when 
-Brot/ll^l — 0.12. Furthermore, Tbar becomes only about 10 
yrs when i?rot/|W^| approaches 0.15. 



Durisen J1977t found that the angular momentum loss time 
scale is comparable to the growth time scale of the insta- 
bility (Tj^bar — Tbar), and therefore we may conclude that 
the bar-mode instability can be an efficient mechanism to re- 
move angular momentum for white dwarfs with £'iot/|W'|> 
(-E'rot/ll^Dbar. crit, as long as this is not suppressed by turbu- 
lent motions (see Sect. 15. 3t . 

5.2. r-mode instability 

Recently, Andersson (1998i and Friedman & Morsink J1998> 
have found that r-modes in all rotating inviscid stars are unsta- 
ble to the CFS instability. The possibility of gravitational wave 
radiation from fast rotating white dwarfs due to the r-mode in- 
stability was also discussed by Hiscock (.1998J . Andersson et 
al. ( fT9^ and Lindblom 



119991 1. Lindblom {199% gives the 
growth time scale for this instability with m = Z = 2 as 
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Fig. 14. Growth time scales for the bar mode instability (rbar, 
Eq. I20> as well as for the r-mode instability (ti, Eq. 12 It . 
Tbar is plotted as function of £'rot/|W^|, with thin lines for 3 
different masses: 1.4 M0 (solid line), 1.5 M0 (dotted line) 
and 1.6 Mq (dashed dotted line). Two different values of 
{Erot/\W\)hai, crit (0.1 and 0.14) for the bar-mode instability 
are considered as indicated in the figure. The thick dashed line 
denotes Ti in the models of sequence A6, while the thick solid 
line is for sequence A 12. The thick dashed-dotted and dotted 
lines give the viscous time scale for sequences A6 and A12 
respectively. 

This time scale is plotted in Fig. for two model sequences 
(A6 and A12). Interestingly, we find that the r-mode instabil- 
ity becomes important when E^ot/\W\<, {Erot/\W\)haT, crit, 
compared to the bar-mode instability. The sensitivity of the 
time scale t^- to £'rot/|W^l is much weaker than that of Tbar, 
and Tr remains in the range 10^ . . . 10^ yr for £'rot/|W^|< 
(-Erot/lM^Dbar, crit- Notc that in our models has much 
smaller values than in Lindblom J1999> . who found to be 
larger than 6 x 10^ yr in the observed DQ Her objects. This 
difference is mainly due to the fact that our white dwarf mod- 
els are rotating differentially, with the outer layer rotating faster 
than the inner core. More precise estimates of Tj require, of 
course, a multidimensional study, which is beyond the scope 
of this paper Nevertheless, the present estimates indicate that 
the r-mode instability might be a promising mechanism for the 
removal of angular momentum in a differentially rotating mas- 
sive white dwarf within an interestingly short time scale (^ 10^ 
yr),whenSrot/|W^| < {E,ot/\W\) 

bar. crit 

5.3. Implications for the final fate of accreting white 
dwarfs 

The previous sections imply that our differentially rotating 
massive white dwarfs can lose angular momentum via the bar- 
mode instability if i5i-ot/|W^| ^ (£'rot/|W^|)bar, crit, and via the 
r-mode instability if Erot/\W\ < (i^rot/| W^Dbar, crit- These 



instabilities would remove angular momentum from the outer 
layers of the white dwarfs — as shown in the numerical simula- 
tions of rotating neutron stars (Lindblom et al. l2002t — where 
most of the white dwarf angular momentum is located in our 
models. 

Although the CFS instabiUty is likely suppressed in the 
presence of strong viscosity for both, the r-mode and the 
bar-mode (e.g. Andersson & Kokkotas 2001 1, Imamura et al. 
J1995t suggest that secularly unstable modes may not be neces- 
sarily damped in the presence of a large effective viscosity due 
to turbulence. Given this uncertainty, we consider both cases, 
i.e., we assume a strong damping of secular instabilities due 
to turbulence induced by the shear instability as Case 1, and a 
persistence of secular modes even when the shear instability is 
present, where the CFS instability is only affected by the mi- 
croscopic viscosity (Case II). 

5.3.1. Case I 

If the unstable CFS modes are damped by turbulence, an ac- 
creting white dwarf will not experience any angular momentum 
loss due to the CFS instability during the mass accretion phase 
in which the dynamical and secular shear instability persist 
throughout the white dwarf interior (see Sect. 14. 2> . no matter 
how large i?i.ot/|M^| becomes. Once the mass accretion ceases, 
the angular momentum redistribution in the white dwarf will 
continue until the degree of the differential rotation becomes 
weak enough for the shear instability to disappear Once the 
fluid motion in the white dwarf becomes laminar, only the elec- 
tron and ion viscosities will contribute to the viscous friction. 

The time scale of viscous dissipation (r^) through the elec- 
tron and ion viscosity is plotted in Fig.[2]for model sequences 
A6 and AlO. Here t^, is calculated following the equation (4.2) 
in Lindblom ( 1999 1. It is found that t^, is far larger than and 

Tbar (except for Thar at E,ot/\W\ ~ (£'rot/|VF|)bar, crit)- This 

implies that the unstable modes of the CFS instability will not 
be damped by the microscopic viscosity before they grow to 
a dynamically meaningful level, once the shear instabihty has 
decayed. 

If the white dwarf mass has grown to A/wd?^ 1-4 M0 by 
the end of the mass accretion phase, carbon ignition at the white 
dwarf center will be delayed as follows. After the mass accre- 
tion stops, the white dwarf will evolve without losing angular 
momentum until the shear instability decays. When the white 
dwarf interior becomes laminar, the white dwarf will begin to 
lose the angular momentum either by the bar-mode instabiUty 
if Erot/\W\ > (i^rot/|T^|)bar, crit, or by the r-mode instabil- 
ity if £'rot/|VF| < (i^rot/|W^|)bar, crit- The corc density wifl 
increase as it loses angular momentum and carbon will ignite 
at the center, eventually. 

Can such a white dwarf end in a Type la supernova? It de- 
pends on how much time is spent from the halt of mass ac- 
cretion to carbon ignition, the time scale to which we will re- 
fer as Tdciay- Nomoto & Kondo il993t found that if the core 
of a white dwarf has solidified, core carbon ignition is likely 
to end in a collapse rather than an explosion. I.e., to obtain 
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Fig. 15. Evolution of the angular velocity in (a), and tempera- 
ture in (b), of a rapidly rotating 1 .50 Mq white dwarf after the 
mass accretion stops. The initial model at t = has been taken 
from model sequence A 12. 



a SN la, carbon ignition should occur before the white dwarf 
core is crystallized. 

In an attempt to estimate rdoiay, we have selected a white 
dwarf model which has grown in mass to 1 .50 M0 in model 
sequence A 12, and let it evolve without further mass accre- 
tion. The evolution of the angular velocity in the white dwarf 
model is shown in Fig.^] The initial model of this calculation 
is characterized by Tc 1.37 x 10*^ K, pc = 1-87 x 10* g/cm^ 
and ii^iot/lM^I = l-H- As the white dwarf cools, the buoyancy 
force in the white dwarf interior becomes weaker and thus the 
dynamical shear instability continues to operate in the central 
core even though the degree of the differential rotation becomes 
continuously smaller. Only when t ~ 1.1 x 10* yr, the shear 
strength becomes so weak that the shear instabihty finally de- 



cays. The central temperature and density at this moment are 
4 x 10'' K and 1.57 x 10* g/cm"^ respectively. This example 
indicates that the CFS instability can begin to operate about 
10* yr after the end of the mass accretion. Which instability 
mode is responsible for the following angular momentum loss 
depends on the value of (fi'rot/lM^Dbai , crit of the white dwarf 
configuration at this point. 

If the angular momentum loss time scale is short enough, 
(i.e., Tj < 10^), the white dwarf may end in a Type la su- 
pernova, since the central temperature is still relatively high, 
and it will increase by the contraction of the core as the white 



dwarf loses angular momentum. Piersanti et al. J2003t showed 
that an initially cool Chandrasekhar mass white dwarf with 
Tc ~ 5.6 x 10^ K is heated to - 10* K if the white dwarf 
loses angular momentum on a time scale of 10^ . . . 10*^ yr 

In the example given above, i?rot/|W^| is decreased from 
0.11 to 0.10 during the evolution. If (i?rot/|W^|)bar, crit at 
this point is less than 0.1, the bar-mode instability will re- 
move angular momentum in Tj bar, which may be short enough 
(i.e., < 10^ yr) to produce a SNe la. On the other hand, if 
(£'rot/|W^|)bar, crit IS larger than 0.1 at the given moment, the 
r-mode instability will become dominant. The value of at 
this point is about 2 x 10^ yr. If the time scale for the following 
angular momentum loss is of the order of Ti , the white dwarf 
will safely end in a SN la. 

Although this discussion is somewhat speculative, it may 
suggest an interesting implication for the masses of white 
dwarfs in SNe la explosions. If the CFS instability is sup- 
pressed by the shear instability, white dwarfs can grow in 
principle until they become dynamically unstable (i?rot/|W^|<^ 

0. 27) without becoming secularly unstable. In most close bi- 
nary systems, the mass accretion will stop before reaching this 
point, given that the mass at Srot/lW^I = 0.18 (A/o.is) in our 
models can be already as high as 1 .45 Mq - 1 .75 Mq (Table|3}. 
Therefore, the upper limit for the white dwarf mass at the mo- 
ment of the supernova explosion (A/sn) is determined by the 
maximum possible mass a white dwarf can achieve by the mass 
accretion, (A/max)- The lower limit for A/sn is simply 1.4 Mq. 

1. e., the white dwarf mass at the moment of supernova explo- 
sion will be inbetween 1 .4 Mq and Afmax (1-4 Mq < AfsN ^ 

Afmax). 

The lower limit for AfgN (1-4 Mq) can not be different for 
different initial masses. However, the upper limit (Afmax) may 
be systematically larger for higher initial masses. I.e., Langer 
et al. ( 2000 1, found that Mmax is larger for a higher A/jnit in 
main sequence star + white dwarf binary systems. Furthermore, 
Afmax is expected to show a dependence on the metalicity of 
the mass donor Langer et al. ( 2000) showed that A/max for 
main sequence star + white dwarf binary systems decreases 
for lower metalicity. This implies a higher probability to ob- 
tain massive exploding white dwarfs for higher metalicity. If 
more massive white dwarfs gave brighter SNe la as Fisher et 



al. J 19991 1 speculate, the brightest SNe la observed at low met- 
alicity may be dimmer than those at higher metalicity, while 
the luminosity of the dimmest ones may not be much different 
(unless affected by the CO ratio; Umeda et al. 1 19991 Hoflich et 
al.. 2000b. 
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5.3.2. Case II 

If the CFS instability is not affected by the shear instability, 
an accreting white dwarf will start losing angular momentum 
when Ti- or Tbar becomes smaller than the accretion time scale 
(race ^ AM/M). 

If we consider sequence A6 as an example, Tj becomes 
comparable to Tacc only when £'iot/|W^| — 0.14. Therefore, 
if we assume (i?iot/|W^|)bar, crit = 0.1, the bar-mode instabil- 
ity will dominate before the r-mode instability becomes im- 
portant. Once £'rot/|VF| becomes larger than 0.1, the white 
dwarf will lose angular momentum via the bar mode instability 
while gaining angular momentum continuously from the ac- 
creted matter. It is likely that Srot/lW^I will not change much 
after Tbar ~ ^acc IS obtained. This means that ii^rot/lW^I can not 
increase much from the critical value of 0.1, since Tbar drops 
rapidly as it deviates from (-Erot/|W^|)bar, crit (Fig. ll4t . If mass 
accretion continues, the white dwarf will finally reach carbon 
ignition at il/wD — 1-76 Mq, according to the results of se- 
quences C9 . . . CI 2, where carbon ignites at the center when 
MwD ^ 1.76 Mq, with E,ot/\W\ ~ 0.11 (see also Sect.|7}. If 
mass accretion stops before carbon ignition is reached but after 
the white dwarf mass has grown beyond 1.4 Mq, carbon igni- 
tion will be delayed until the white dwarf loses enough angular 
momentum. 

In summary, we can expect that even in Case II the mass of 
exploding white dwarfs will show a significant diversity, as in 
Case I. The lower limit for Msn should be ^^1.4 Mq in both 
cases, but the upper limit may be either the critical mass for 
carbon ignition at that E'rot/jW^I which gives tj « Tacc, or the 
maximum possible achievable mass by mass accretion. 

5.3.3. Other uncertainties 

In reality, the mass transfer rate in a close binary system does 
not remain constant. For instance, in a binary system consist- 
ing of a main sequence star + a white dwarf, where thermally 
unstable mass transfer occurs, the mass transfer rate increases 
rapidly from the onset of the mass transfer, followed by slow 
decrease to the point that the nuclear shell burning can not be 
stable any more (Langer et al. 2000). Although shell sources 
may be more stable with rotation (Yoon et al. l2004i) . the white 
dwarf will experience strong shell flashes if the mass accretion 
rate decreases below 10^^ Mo/yr. This may cause a signif- 
icant loss of mass, which may even decrease the white dwarf 
mass. This may also lead to a removal of angular momentum 
from the white dwarf (Livio & Pringle [T998l l. Effects of mag- 
netic fields, if they are strongly amplified due to the differen- 
tial rotation (Spruit 2002 1, may also serve to break the white 
dwarf by magnetic torques and/or by magnetic dipole radia- 
tion. These possibilities will be subject of further studies. 

6. Implications for the SNe la diversity 

We note that for all considered circumstances, effects of rota- 
tion in accreting white dwarfs offer the possibility of SN la 
explosions from white dwarfs more massive than ^1.4 Mq 
(hereafter designated as super-Chandrasekhar mass). The ob- 



servational signatures of super-Chandrasekhar mass explosions 
are currently difficult to predict. The peak brightness and light 
curves of SNe la might be a sensitive function of the igni- 
tion conditions such as the CO ratio, the density at the mo- 
ment of carbon ignition, the speed of the deflagration front, the 
core rotation rate, and possibly the transition density at which 
the deflagration turns into a detonation (e.g. Khokholov l 1 99 1 al 
1991bi Hoflich & Khokhlov 1996 Niemever & Wooslev [T997l 
Hoflich et al. 1998 Iwamoto et al. 119991 U meda e t al. [T999l 
Hillebrandt & Niemeyer 2000 Woosley et al. l2003l . 

If the ignition conditions are not strongly affected by the 
mass of the exploding white dwarf, it is likely that more mas- 
sive one gives a brighter SN la, since a more massive white 
dwarf can provide more fuel to produce ^^Ni. For instance, pro- 
genitors with a super-Chandrasekhar mass are often invoked to 
explain such an anomalous SN la as SN 1991T. Fisher et al. 
J 19991 1 note that the mass of ^®Ni produced from the explosion 
models of a Chandrasekhar mass white dwarf remains to be 
smaller than about 0.9 Mq even in the case where a pure deto- 
nation is considered (e.g. Hoflich & Khokhlov' 1996V Although 
the nickel masses derived from the luminosities of normal SNe 
la are constrained to 0.4 .. . 0.8 Mq (Leibundgut 2000), the 
existing models fail to explain such a peculiar SN la as 1991T, 
whose brightness implies the production of about 1.0 Mq 
of ^^Ni. The fact that explosion models with the canonical 
Chandrasekhar mass could not explain the peculiarity of SN 
199 IT suggests that it is worthwhile to investigate the possible 
outcome of super-Chandrasekhar mass explosions, even if we 
can not exclude the possibility that the luminosity of 199 IT is 
biased by the uncertainty in determining the distance of its host 
galaxy (Hanato et al. l2002t . 

On the other hand, Fisher et al. J1999> suggested the ex- 
plosion of a super-Chandrasekhar mass white dwarf from the 
merger of double CO white dwarfs as explanation of the over- 
luminosity found in SN 199 IT. We note that differentially ro- 
tating single degenerate progenitors may be a more natural ex- 
planation for the super-Chandrasekhar mass scenario than dou- 
ble degenerate mergers, which fail to produce SNe la in nu- 
merical models due to the off-center carbon ignition induced 
by the fast mass accretion with M > 10^^ Mq/jt (Saio & 
Nomoto [T985l[T998t . 

Rotation also bears implications for the polarization of 
SNe la. According to our results, the rotation velocity of ex- 
ploding white dwarfs may strongly depend on the history of 
angular momentum loss via gravitational wave radiation. The 
more angular momentum is lost, the slower the white dwarf 
rotation will be, which means that the polarization strength in 
SNe la may show a significant diversity. Some white dwarfs 
may have a chance to end in an SN la while they are ro- 
tating very rapidly, as in sequences C9, CIO, Cll and CI 2, 
where carbon ignition occurs at the white dwarf center when 
AfwD ^ 1-76 Mq andi^rot/lW^I ~0.11.Such exploding white 
dwarfs will show strong features of asphericity in the explo- 
sion, which might give a plausible explanation for the polariza- 
tion observations in SN 1999by (Howell et al. l2001> and SN 
2001 el (Wang et al. .2003> Kasen et al. .20031 . 
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7. Pre-explosion conditions of a fast rotating 
whiite dwarf 

Rotation may have interesting consequences for the supernova 
explosion itself. Fig. ^] shows the physical properties of the 
white dwarf model of sequence C9 when AfwD = 1-76 Mq, 
where the central temperature and density reach 5.5 x 10* K 
and 2.2 x 10^ g/cui^, respectively. The central region (0.0 ^ 
M). ^ 0.6Mq) is found convectively unstable due to carbon 
burning, rotating rigidly due to the convective angular mo- 
mentum redistribution. The thermonuclear runaway is expected 
to develop when the central temperature reaches a few times 
10^ K. The convective core will be more extended by then, 
spinning-up the central region further In the model shown in 
Fig. the central region rotates with w^ot — 1000 km/s, 
which is well above the convective velocity (~ 100 km/s; e.g. 
Hoflich & Stein 2002i and somewhat larger than the expected 
initial velocity of the deflagration front (wdcf .init — O.OSwsound, 
e.g. Nomoto et al. 1984 1. The equatorial rotation velocity in- 
creases sharply from the center to the edge of the convective 
layer (M, = 0.6 M©), by a factor of 10. 

This fast rotation in the central carbon burning core may 
affect the explosion as follows. Numerous studies indicate that 
the combustion front develops large scale turbulent motions, 
especially due to the Rayleigh-Taylor instability (e.g. Reinecke 
et al. 120021 Gamezo et al. 120031 and references therein). It is 
a crucial question whether this turbulence can induce a deto- 
nation (Khokhlo v J991a J991b, Niemeyer & Woosley 19971 
Khokhlov et al. [T9971 . Alfliough pure deflagration models 
show a good agreement with observations (Nomoto et al. l 19841 
Hoflich & Khokhlov 1996j, recent three dimensional calcula- 
tions show that large amounts of unburned and partially burnt 
materials are left even at the center (Reinecke et al. 120021 
Gamezo et al. 2003 1, which is not supported by observations 
(however, see Baron et al. I2003> . This problem would dis- 
appear if a detonation were triggered by the turbulent defla- 
gration ('delayed detonation', Khokhlov 199 la i. Furthermore, 
delayed detonation models are shown to reproduce the light 
curves and spectra of well observed Type la SNe (Hoflich & 
Khokhlov 1996 Hoflich et al. 1996 1. A study of the nucleosyn- 
thesis output of SNe la by Iwamoto et al. (119991 1 also favors the 
delayed detonation scenario. 

Several authors conclude that the probability for triggering 
a detonation can be enhanced if the turbulence is strong enough 
to form a thick mixed region of unburned and burning mate- 
rials with a constant temperature gradient (Khokhlov ^199 lb 



Niemeyer & Woosley 119971 Khokhlov et al. 119971 Lisewski 
However, no robust mechanism for the forma- 



et al. I200m 

tion of such a strong turbulence has yet been suggested. 
(Niemever [T999t . 

A rapidly rotating white dwarf like the one shown in Fig. 1161 
will be spun down due to expansion once the thermal runaway 
occurs, as the deflagration wave propagates. However, the outer 
layers will still rotate rapidly at Vrot ^ 10^ km/s. Lisenwski et 
al. (2Q00l show that strong turbulence, i.e. wturb ~ 10'^ km/s, is 
necessary for a detonation to form successfully. Since the rota- 
tion velocity in the outer region will be of this order, the shear 
motion may provide enough kinetic energy for the turbulence 
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Fig. 16. (a) Temperature (solid line) and angular velocity 
(dashed line) as function of the mass coordinate in the last com- 
puted model of sequence C9. (b) Equatorial rotational velocity 
(solid line) for the same model as in (a) as function of the mass 
coordinate. The dashed line denotes the sound speed as func- 
tion of the mass coordinate. 



intensity to satisfy the condition for a transition of the defla- 
gration into a detonation, which might not be possible in the 
non-rotating case since the maximum value for Wturb is ~ 10^ 
km/s in the Rayleigh-Taylor-driven turbulence. Alternatively, 
the fast rotation in the outer region can tear the turbulent defla- 
gration front and enhance the nuclear burning surface, which 
may be another way to trigger a detonation (Hoflich 2003, 
private communication). Future multi-dimensional calculations 
are required to test the validity of these scenarios. 



8. Gravitationai wave radiation in pre-explosion 
whiite dwarfs? 

Rotation may have another important and exciting observa- 
tional consequence. Our results indicate that the white dwarfs 
in SN la progenitor systems may emit gravitational wave radi- 
ation (GWR). Therefore, if we could observe GWR preceding 
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Table 4. Expected properties of the gravitational waves from 
the selected accreting white dwarfs. 



No. 


M 


E,at/\W\ 


/ibar at 10 MpC 


/gw 




Mq 




10-24 


Hz 


A2 


1.18 


0.10 


1.7 


0.22 


A6 


1.30 


0.10 


2.9 


0.32 


AlO 


1.42 


0.10 


4.3 


0.46 


A2 


1.34 


0.14 


3.7 


0.32 


A6 


1.44 


0.14 


5.0 


0.43 


AlO 


1.55 


0.14 


7.5 


0.62 



a SN la event, this would provide strong evidence for the sce- 
narios outlined in Sect.|5] 

The GWR due to the CFS bar-mode instability may be 
much stronger than the r-mode signal, since perturbations due 
to the bar-mode involve larger density changes than those due 
to the r-mode (e.g. Andersson & Kokkotas 2001). Fryer et 



al. (I2002> derive the gravitational wave amplitude observed at 
Earth from a source which undergoes the bar-mode instability 
at a distance d as 



1/2 



G MR^VL^ 



(22) 



where R denotes the equatorial radius of the source, M the 
mass and Q, the angular velocity. The frequency of the gravita- 
tional waves (/gw) due to the bar mode instability is twice the 
rotation frequency of the source. We estimate /ibar and /gw for 
selected white dwarf models in Table@] For these calculations, 
we used mean values for i7 and R given in Table |2land|3] 

Table |4] indicates that /gw is within the range of 
0.1 . . . 1.0 Hz. Observation of the GWs with such frequencies 
may be performed by low frequency gravitational wave detec- 
tors such as LISA, which will cover the range 10-^ — 1.0 Hz 
(e.g. Cutler & Thorne "2002'l. The expected strength of the 
gravitational wave signal is about IQ-^^ ai d = 10 Mpc. 



Hiscock (I1998> estimated the strength of GWs due to the r- 
mode instability from white dwarfs in the observed DQ Her 
systems as h ~ lO-^'^. His calculations show that this strength 
is well above the detection limit of the LISA interferometer. 
Therefore, the results shown in Table |3 imply that SNe la 
progenitors in nearby galaxies could be within the observable 
range. Given that GWs from SN la progenitors will be emitted 
over a secular time, and that SNe la are observed with a rela- 
tively high rate of 3 x IQ-'^ yr^^ per galaxy, the probability 
to detect the GW signal may be significant. 

Even if a white dwarf would not reach the bar-mode in- 
stability at (i?rot/|W^|)bar, crit, it might Still emit gravitational 
waves via the r-mode instability. Although their strength is 
rather unclear and may be much weaker than that of the bar- 
mode, they may still be detectable if the source is close enough. 

9. Conclusion 

We summarize the results of this paper as follows. 

1 . The role of the Eddington sweet circulation, the GSF in- 
stability and the shear instabiUty for the transport of angular 



momentum in non-magnetized white dwarfs has been investi- 
gated (Sect.|2j- Although Eddington sweet circulation and the 
GSF instability are important for the redistribution of angular 
momentum in the non-degenerate envelope, their importance is 
small in the degenerate core compared to the shear instability. 
The secular shear instability can not operate in the strongly de- 
generate core, since the thermal diffusion time becomes longer 
than the turbulent viscous time for densities higher than a crit- 
ical value (i.e., p > 10® . . . 10^ g/cm^, Fig.|2li. On the other 
hand, the criterion for the dynamical shear instability is signif- 
icantly relaxed for higher density because the buoyancy force 
becomes weaker with stronger degeneracy (Fig.[0. As a result, 
the degenerate inner core is dominated by the dynamical shear 
instability in accreting white dwarfs as shown in Sect. 14.21 

2. We have followed the redistribution of the angular mo- 
mentum in accreting white dwarfs by the above mentioned pro- 
cesses ("Sect. l4~2l . We find that accreting white dwarfs do not 
rotate rigidly, but differentially throughout their evolution, for 
the considered accretion rates {M = 3 ... 10 x 10-^ Mo/yr). 
In the degenerate core, once the shear factor decreased to the 
threshold value for the onset of the dynamical shear instabil- 
ity, the time scale for further angular momentum transport be- 
comes larger than the accretion time scale. Accordingly, strong 
differential rotation is retained in the inner core with a shear 
strength near the threshold value for the dynamical shear insta- 
biUty (Fig.|7}. 

3. Accreting white dwarfs, as they rotate differentially, may 
not reach central carbon ignition even when they grow be- 
yond the canonical Chandrasekhar limit of ~1.4 M©. This is 
in accordance with previously obtained results based on other 
methods (Ostriker & Bodenheimer fTgeSlfTW ? Durisen [T975bl 
Durisen & Imamura 198 1 1. A secular instability to gravitational 
wave radiation through the r-mode or the bar-mode may be im- 
portant for determining the final fate of accreting white dwarfs. 
The masses of exploding white dwarfs are expected to vary in 
the range from the canonical Chandrasekhar mass (~1.4 M©) 
to the maximum possible mass that the white dwarf can achieve 
by mass accretion in a binary system (~ 2.0 M©, Langer et 
al. l2000> . This may have consequences for the diversity in the 
brightness and polarization of SNe la. 

4. Fast rotation in the white dwarf core may change the 
supernova explosion since it can affect the evolution of the tur- 
bulent nuclear burning flames by providing a large amount of 
turbulent kinetic energy and/or by enhancing the burning sur- 
face significantly. It needs to be clarified whether this may be a 
plausible mechanism to induce the transition from deflagration 
to detonation. 

5. White dwarfs which accreted enough angular momen- 
tum may be detectable sources of gravitational waves in the 
near future. Our models show that these will emit gravitational 
waves with frequencies of 0.1 — 1.0 Hz. Space-based interfero- 
metric gravitational wave detectors such as LISA could observe 
such signals from rapidly rotating SNe la progenitors in nearby 
galaxies. 
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